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Abstract 

We show that the twist-2 light-cone wave- functions of r/ c and J/ip can be factorizcd with nonrela- 
tivistic QCD (NRQCD) at one-loop level, where the nonperturbative effects are represented by NRQCD 
matrix elements. The factorization is achieved by expanding the small velocity v, which the c- or c- 
quark moves inside a rest quarkonium with. At leading order of v the twist-2 light-cone wave-functions 
of rj c and J/ ip can be factorized as the product of a perturbative function and a NRQCD matrix ele- 
ment. The perturbative function is calculated at one-loop level and free from any soft divergence. Our 
results can be used for the production of J/i/jr] c through e + e~ -annihilation and of a charmonium in 
B-decays, which are studied in experiment of two £?-factories. 
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Exclusive processes involving a charmonium can be studied with nonrelativistic QCD(NRQCD)pQ. 
One uses the fact that the c- or c-quark moves with a small velocity inside a rest charmonium and makes 
an expansion in v to factorize the perturbative- and nonperturbative effects. In this approach of NRQCD 
factorization the production of a charmonium can be viewed as a two-step process, where a free cc-pair 
is first produced and then the pair is transformed into the bound state of the charmonium. Because 
the charm quark mass m c is large, the production of the cc-pair can be calculated with perturbative 
QCD. The transformation is nonperturbative and can be characterized with NRQCD matrix elements. 
Applications of the approach can be found in the comprehensive reveiwj^j. 

Although the approach of NRQCD factorization is consistent in theory, but in practices there are 
problems when confronting experiments. First, the expansion in v converge slowly because v for char- 
monia is not small enough with the typical value v 2 ~ 0.3. Secondly, since in the approach the charm 
quark mass must be kept as nonzero, there will be in the perturbative expansion large logarithms terms 
like ln(m c /Q), when the process of the production involves some large scales, generically denoted as Q, 
and Q >> m c . These large logarithms terms will spoil the applicability of perturbative QCD. 

Alternatively, one can use the collinear factorization for those exclusive processes involving some 
large scales with Q >> m c . It should be noted that the collinear factorization for exclusive processes of 
light hadrons has been suggested long time ago[31|3]. In this approach the nonperturbative effects are 
parameterized with different light-cone wave-functions(LCWF's) of a charmonium and the amplitude of 
a process can be written as a convolution with those LCWF's and a perturbative function. In calculating 
the perturbative function charm quarks should be taken as light quarks at leading twist, i.e., m c = 0. 
Therefore, large logarithms terms ln(m c /Q) will not appear in the perturbative function, but will appear 
in LCWF's. Then those large logarithms terms can be easily re-summed by using renormalization group 
evolutions of LCWF's. The situation here is similar to inclusive production of a charmonium, where one 
uses the concept of fragmentation functions to resum large logarithms terms. 

Experimentally exclusive production of a charmonium in S-decays have been studied extensively 
at two £?-factories. An interesting study at U-factories reveals a puzzle for theory. The cross-section 
of e + e~ — > J/ipr] c has been measured [5]. The measured cross-section is about an order of magnitude 
larger than theoretical predictions. The theoretical predictions are based more or less on the approach 
of NRQCD factorization |0]. Adding one-loop correction in the approach[Jj, the discrepancy still remains 
although the correction is significant. It has been suggested [51 1§1 ITU] that by using the collinear factoriza- 
tion it is possible to solve the puzzle. In this approach the nonperturbative effects related to charmonia 
are parameterized with different LCWF's of charmonia. The collinear factorization has been also used 
in other exclusive processes involving quarkonia|llj. 

In this letter we point out that LCWF's of charmonia can be studied with NRQCD, in which one can 
factorize perturbative- and nonperturbative effects. The perturbative effects are related to the fact that 
charm quarks are heavy. The situation is similar to parton fragmentation functions into a charmoium 
for inclusive production, where a NRQCD factorization can be performed [12]. We will show here that 
the NRQCD factorization for LCWF's can be performed at one-loop level for twist- 2 LCWF's of J/tp or 
rj c , i.e., of a S- wave quarkonium, where we make an expansion in v and only keep the leading order. At 
the order the shape of LCWF's is completely determined by perturbative QCD. It should be noted that 
LCWF's of charmonia have been studied in a potential model[T5]. 

We will use the light-cone coordinate system, in which a vector is expressed as = (a + , a~, aj_) = 
((a + a 3 )/V2, (a - a 3 )/V2, a 1 , a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . There is one LCWF of r? c at leading twist, 
which is defined as: 



</> Vc (u, (J,) 



I 



dx 



e+iu p+x ( |c(0)Lt (oo, 0) 7 +75^n(oo, x-n)c(x-n)\ Vc (P)) 



(1) 
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where ry c carries the momentum = (P + ,P ,0,0). L n is the gauge link along the light-cone direction 
= (0, 1, 0, 0) to make the definition gauge- invariant: 



L n (oo,z) 



Pexp(-ig s J dXn ■ G(Xn + z)\ . (2) 



The variable u is the momentum fraction carried by the c-quark, i.e., the c-quark has the +- component 
uP + of its momentum. Similarly, the LCWF's of J/tp with different polarizations A are defined as: 

0[| = / ^-e +iuP+x ~ (0\c(0)Ll(oo,0h + L n (oo,x-n)c(x'n)\J/ij(P,X = 0)), 

<f>±(u,») = |^e +iuP+ ^(0|c(0)LUoo,0)(-ia + ^(A ± ))L n (oo,x-n)c(x-n)|J/^(P,A ± )), (3) 

where is the polarization vector of J/tp andn-e*(A = ±1) = l-e*^(X = ±1) = 0. (e* (A = ±1) = e*(Aj_). 
There are only two LCWF's of J /v/j at leading twist. These definitions can be easily obtained from those 
of light hadrons |14M15| . the definitions here are only different in the normalization than those defined in 
|14l I15j . The LCWF's are Lorentz-invariant along the z-direction. 

In the framework of NRQCD based on the small velocity v expansion, one expects that those functions 
take a factorized form: 

— \ / 1 1 hp * ' 



where 4> n {n = rj c , ||, _L) can be calculated in perturbative QCD and they are free from any infrared- and 
Coulomb singularity, as we will show. They are dimensionless. H stands for J/i/j or r/ c . The matrix 
element (0\Oh\H) is defined with NRQCD operators, whose definitions will be given later. They contain 
all nonperturbative effects. The factor appears because the states in the definitions of LCWF's 

are covariantly normalized, while the states in NRQCD matrix elements are usually nonrelativistically 
normalized. 

To extract the perturbative function <p n , we replace the quarkonium state with a cc-pair state. After 
the replacement, one can calculate LCWF's in the left hand side of Eq.(4) and the matrix elements in the 
right hand side with perturbative QCD and then extract <p n . In such a calculation, infrared- and Coulomb 
singularities can appear in the left hand side, but the same will also appear in the matrix elements in the 
right hand side so that <j) n is free from any soft divergence. This will be shown in our calculation below. 
Because LCWF's are Lorentz invariant in the ^-direction, one can calculate LCWF's in the rest frame of 
the pair. In the rest frame the cc-pair state is specified with \c(p±)c(p2)} , the momenta are given by: 

pi = m c v, p 2 = -m c v, P = Pl +p 2 (5) 

For easy calculation we let the quarks move in the z-direction, i.e., v = (0, 0, v) and v > 0. This will not 
bring any problem for S- wave quarkonia at leading order of v. For convenience we introduce the variable 
uj the following: 

" = |b u>= l -+ l -v + 0{v% (6) 



and denote 

r — r 

We will use Feynman gauge in this work. 



{7 + 75, 7 + , -™ + %}- (7) 
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We will always use the notation for a quantity A as A = A^ ' + Ay^ + • ■ ■, where A^ denotes the 
tree-level contribution and A^ denotes the ra-loop contribution. At tree-level we have: 

<t>£\u) = 8{u-u)-j^v{p 2 )T nU { Pl ). (8) 

The product of the spinors in the above can be expressed with two-component NRQCD spinors £ and 77 
in the limit v — > 0: 

v{P2h + 75u{ P i) = V2m c rft + 0(v 2 ) = V2m c (0\ X ^\c(pi)c(p 2 )) {0) + 0(v 2 ), 
v(P2h + u( Pl ) = V2m c rfa 3 £ + 0(v 2 ) = V2m c {0\x^ 3 i>\c( Pl )c(p2)) {0) + 0(v 2 ), , 
v(P2) (-^elj u(pi) = V2m c rfa ■ e*^ + 0(v 2 ) = y/2m c (Q\x}<r ■ e* ± ^\c( Pl )c(p 2 )) i0) + 0(v 2 ), (9) 

where ip and % are NRQCD field operators, ip annihilates a c-quark and x creates a c-quark. With the 
above expression we can define two NRQCD matrix elements relevant to our work: 

(o\o Vc \ Vc ) = (o\xH\vc), (o\o m \JM = <o| x tI]r<x-e*(A)|JMA)>. ( 10 ) 
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In the above the quarkonium is in the rest frame. With the above results we can specify the operators 
in Eq.(4) with <j) n determined at tree- level: 

0\\,±( u > N = 0||,±(«>aO T~F= + 

= ^\u^) = ^\u^) = 5{u-l/2). (11) 



Figure 1: The virtual part of the correction. The double lines are for the gauge link. 

An important purpose of our work is to show that the perturbative functions can be calculated with 
perturbative QCD at one-loop level and they are free from any soft divergence. For this purpose we need 
to calculate 4> n and the NRQCD matrix elements at one-loop to extract 4> n , where we take the partonic 
state \c(pi)c(p2)) ■ As usual, the one-loop correction to 4> n can be divided into a virtual- and a real part. 
The virtual part is represented by diagrams given in Fig.l. The contributions from Fig. la and Fig. lb are 
from the external legs. They read: 

<Pn\la(u) + <t>n\lb{u) = Z w $Xu), (12) 

with 



2as 

3ir \ ei fij ' 2 fi' 2 



Z w = =^(-+ 7 -ln4^ + ln^ + iln^-2 , (13) 
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with ej = d— 4 for the infrared singularity and the scale \i\ for the singularity. The U.V. poles in e = 4 — <i 
are subtracted. The contribution from Fig.lc and Fig. Id reads: 



0n|ic(«) 
0n|ld(«) 



-^#(«)r(|) 



/4V N 


i" 


" 1 l " 


V m c , 




[ej 1-eJ 




r 


" 1 1 " 


V m c - 




Ui 1-eJ 



(14) 



Where T(e/2) is for U.V. divergence and needs to be subtracted. In the above there are also I.R. 
singularities. They will be canceled by those in the real correction. 





Figure 2: The real part of the correction. The double lines are for the gauge link. 



The real correction comes from diagrams given in Fig. 2. The contribution from Fig. 2b and Fig. 2c 
reads 



2\ f / 2 r 



4>n\2b(u) 
<t>n\2c{u) 



^v{ P2 )T n u{ Pl )T{^) J 
2 f A 2\ e / 2 



2(u 



1 



1 



-1+61 



9{u- 1/2) +0(v 2 ), 



9(1/2 -u) + 0(v 2 ). 



(15) 



In the above contributions there are poles at u = 1/2 when e/ = 0. These poles are infrared one and 
should be regularized as the poles in ej. Using the expansion for (^ — u)~ 1+e ' and (u — \)~ l+ei as 
distributions: 



{n -\r^-e ( u- m = I(i)% ( „- 1/2) + ()( „- 1/2 )((^^|i 



(I_„)-i+.,«(l/2-„) 



with the +-distribution defined as: 



u- 1/2 



; ^»/w(^| 

ci V 1/2 - n 



u 



>(rt.)-/(W)(^ 
1/2 dll(/(u) _ /(1/2)) ((lZ^ 

ci V 1/2 - u 



(16) 



(17) 
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where f(u) is a test function. With these results one can show that the the infrared pole in ej is canceled 
in the sum <f) n \i c + <f> n \2c- The same also happens for the sum (f> n \id + (fi n \2b- After subtracting the U.V. 
poles we have: 



4>n\ld{u) + 4>n\2b(u) 



<f>n\lc(u) + <p n \2c(u) 



-2S(u - 1/2) - 40(2u - 1) ln(2u-l 

6 n |i d (l - U) + 4>n\2b{l ~ U) 



5{u - 1/2) -6{2u- I 

ln(2n - 1 
2u - 1 



2«- 1 



(18) 



The most difficult part is the contribution from Fig. 2a. It contains not only U.V.- and I.R. poles, 
but also the Coulomb singularity when v goes to zero, represented as 1/v. After a tedious calculation we 
have: 



'n|2a 



(«) = 



2a s 
3vrP+ 



^(p2)r„u(pi 



2 



2I(u) 



+ 0n|2a,_R(-u), 

(19) 



where (j) n \2a,R,(u) contains a U.V. pole and its form depends on the structure of T n . The contributions 
with the I.R. and Coulomb singularity are universal, they do not depend on T n . The distribution I(u) is 
defined as: 



i(u) 

f 1 9(1 - 2u) 



l-2u) 9(1 -2u) 0(2u-l) 6(2u-l) 



+ 



(l-2u) + (l-2u) 2 + (2u-l) + (2u-l) 
1/2 iu /W-/ (V2) + (V2 - n)/'(l/2) 



2 ' 



JO 



du 



0(2u-l) 
(2« - 1) 2 + 



/(«) 



1/2 



(1 - 2u) 2 
/( U )-/(l/2)-( U -l/2)/'(l/2) 
(2« - l) 2 



(20) 



where f(u) is a test function and f'(u) = df(u)/du. The remainder part is: 



»r) c \2a,R 



(«) = 



2a 



4a,. 



3vrP+ 
+0(2u-l)(l-u) 

4>±]2a,R{u) = 0. 



v(p2)j u(pi) < 9(1 — 2u)u 



lnAr-3-21n(2u-l) 



In 



/i 



-3-21n(l-2u) 



(21) 



To verify the factorization in Eq.(ll) one needs to calculate one-loop correction to the NRQCD 
matrix elements. The one-loop corrections come from the diagrams given in Fig.3. The calculation is 
straightforward. We have: 



<0|O H |c(v)c(-v)> 



(i) 



2a s 
3vr 



'2 , 4ttjU?\ tt 2 



<O|0 H |c(v)c(-v)> 



+ 



2a, (2_ 
3vr I ej 



+ 7 



-\n^-) (0|O^|c(v)c(-v))(°) + O(t; 2 ), 



fj, 



(22) 
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Figure 3: One-loop correction for NRQCD matrix elements. 



The first line is for the contribution Fig. 3a, the second line is for the sum from Fig. 3b and Fig. 3c. With 
above results we can now perform NRQCD factorization and see how those I.R. singularities and Coulomb 
singularity are canceled. 

Comparing the one- loop results of NRQCD matrix elements with the one- loop corrections to cp n , we 
can see that the Coulomb singularity and I.R. singularity from Fig. 2a in <fi n are canceled by the same 
singularities from Fig. 3a in NRQCD matrix elements. The I.R. singularities from Fig. la and Fig. lb in 
4> n are canceled by those from Fig. 3b and Fig. 3c in NRQCD matrix elements. It should be noted that 
the cancelations happen in the diagram-by-diagram way. This leads to expect that the factorization can 
hold beyond one-loop level. After these cancelations, the extracted perturbative functions will be free 
from any soft divergence. We have: 

,2 




(1 - 2u). 

These are our main results. 

From our results we can also derive the dependence on the renormailization scale [i. For T r , 
and 7 + 75: 



(23) 



din fi 2 



2a. 
2vr 



^(1/2 



u) + 20(1 - 2u) 



1 



+ u - 1 + (u 



1 



(24) 



(1 - 2u). 

These LCWF's satisfy the Efremov-Radyushkin-Brodsky-Lepage evolution equation |16j. The equation 
reads: 



d In /x 2 



dvV{u,v)4> nci \\(v,n) 



(25) 
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with 

V(u,v) = 
The +-prescription is 



2a s 
3vr 



1-u ( 1 
- 1 + 

1 — V \ u — V 



u ( 1 

(«-«) + - 1 + 

V V f — M 



9(v - u) 



[V(u,v)]+ = V(u,v) -5(u-v) / dtV(t,v). 

Jo 



(26) 



(27) 



Our results agree with the RG equation. Form our result it is easy to derive the renormailzation group 
equation for cj)±. We have 



d(p±(u,n) 
din fi 2 

V±(u,v) 



dvV±(u, v)cj)±(v, /j.) 
u 



o 

V{u,v) 
3vr 



2a s 
3vr 



0(v -u) + 



1 



u 



1 



[U 



-5(u - v) +B{v 



V 

1 



v) 



{v — u). 



1 



(28) 



where the +-prescription acts on the variable u. 

To minimize possibly large perturbative corrections our results can be used as initial conditions at 
[i = m c to obtain LCWF's through renormalization group equations at a scale \i relevant in a process. 
Through this step the mentioned large logarithms terms ln(Q/m c ), mentioned at the beginning, are 
resummed. In this work, the shape of LCWF's is completely determined by perturbative QCD. Since 
v 2 is not small enough for charmonia, nonperturbative effects can be at higher orders of v 2 can be 
important. However, it is possible to resum these effects. E.g., one can consider exchanges of multiple 
Coulomb glouns in Fig. 2a. The resummation of these exchanges leads to solve a Schrodinger equation. 
The resummation results in that the sharp distribution 5{u — 1/2) becomes smeared. This can be studied 
with potential models as in |13j . 

To summarize: We have shown that in the framework of NRQCD the twist-2 LCWF's of a 5-wave 
quarkonium can be factorized as a product of perturbative functions and a NRQCD matrix element in the 
nonrelativistic limit. The nonperturbative effect is contained in the matrix element and the perturbative 
functions can be calculated with perturbative QCD. They are determined at one-loop level and free 
from I.R.- and Coulomb singularity. The LCWF's have a wide application in exclusive production of 
quarkonia. Our results can be used as initial conditions at [x = m c to obtain the LCWF's at other scale [i 
relevant in a process. With the study performed for twist-2 LCWF's we expect that such a factorization 
also holds for twist-3 LCWF's. A detailed study of twist-3 LCWF's and phenomenological applications 
of these results are under way. 
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